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Abstract: We study the threshold conditions of spatial self organization 
combined with collective coherent optical backscattering of a thermal 
gaseous beam moving in a high Q ring cavity with counter propagating 
pump. We restrict ourselves to the limit of large detuning between the 
particles optical resonances and the light field, where spontaneous emission 
is negligible and the particles can be treated as polarizable point masses. 
Using a linear stability analysis in the accelerated rest frame of the particles 
we derive an analytic bounds for the selforganization as a function of 
particle number, average velocity, temperature and resonator parameters. 
We check our results by a numerical iteration procedure as well as by direct 
simulations of the N-particle dynamics. Due to momentum conservation the 
backscattered intensity determines the average force on the cloud, which 
gives the conditions for stopping and cooling a fast molecular beam. 

© 2008 Optical Society of America 
OCIS codes: (000.0000) General. 



References and links 

1. S. Chu, "Nobel Lecture: The manipulation of neutral particles," C. Cohen-Tannoudji, "Nobel Lecture: Manipu- 
lating atoms with photons," and W. D. Phillips, "Nobel Lecture: Laser cooling and trapping of neutral atoms," 
Rev. Mod. Phys. 70, 685-741 (1998). 

2. E.A. Cornell and C. E. Wieman, "Nobel Lecture: Bose-Einstein condensation in a dilute gas, the first 70 years 
and some recent experiments ," Rev. Mod. Phys. 74, 875-893 (2002). 

3. P. Domokos and H. Ritsch, "Mechanical effects of light in optical resonators ." J. Opt. Soc. Am. B 20,1098-1 130 
(2003). 

4. A. Beige, P. L. Knight, and G. VitieUo, "CooUng many particles at once," New J. Phys. 7, 96 (2005), 
Ihttp: / /www. iop.org/EJ/art icle/ 1367 -2 630/7 /I/ O 96/njp5- 1 -0 96. htmll 

5. V. Vuletic and S. Chu, "Laser cooling of atoms, ions, or molecules by coherent scattering," Phys. Rev. Lett. 84, 
3787-3790 (2000). 

6. H. W. Chan, A. T. Black, and V. Vuletic, "Observation of collective-emission-induced cooling of atoms in an 
optical cavity," Phys. Rev. Lett. 90, 063003 (2003). 



7. P. Maunz, T. Puppe, I. Schuster, N. Syassen, P. W. H. Pinkse.and G. Rempe, "Cavity cooling of a single atom," 
Nature London 428, 50-52 (2004). 

8. S. Nussmann, K. Mun; M. Hijlkema, B. Weber, A. Kuhn, and G. Rempe, "Vacuum-stimulated cooling of single 
atoms in three dimensions," e-print quant-ph/0506067 

9. A. T. Black, H. W. Chan, and V. Vuletic, "Observation of collective friction forces due to spatial self-organization 
of atoms: from rayleigh to bragg scattering," Phys. Rev. Lett. 91, 203001 (2003). 

10. J. Klinner, M. Lindholdt, B. Nagorny, and A. Hemmerich, "Normal mode splitting and mechanical effects of an 
optical lattice in a ring cavity," Phys. Rev Lett. 96, 023002 (2006). 

11. B. Nagorny, T. Elsasser, H. Richter, A. Hemmerich, D. Kruse, C. Zimmermann, and P. Courteille, "Optical lattice 
in a high-finesse ring resonator," Phys. Rev. A 67, 031401(R) (2003). 

12. H.L. Bethlem, G. Berden, F.M.H. Crompvoets, R.T. Jongma, A. J. A. van Roij, and G. Meijer, "Electrostatic 
trapping of ammonia molecules," Nature London 406 491^94 (2000). 

13. N. Vanhaecke, W.D. Melo, B.L. Tolra, D. Comparat, and P. Pillet, "Accumulation of cold cesium molecules via 
photoassociation in a mixed atomic and molecular trap," Phys. Rev. Lett. 89, 063001 (2002). 

14. R. Bonifacio, C. Pellegrini, and L.M. Narducci, "Collective instabilities and high-gain regime in a free electron 
laser," Opt. Comm. 50, 373-378 (1984). 

15. R. Bonifacio, 1. De Salvo, L.M. Narducci, and E.J. D' Angelo, "Exponential gain and self-bunching in a collective 
atomic recoil laser," Phys. Rev. A 50, 1716-1724 (1994). 

16. J. Guo, P. R. Berman, and B. Dubetsky, "Recoil-induced resonances in nonhnear spectroscopy," Phys. Rev. 
A 46,1426-1437 (1992) and "Comparison of recoil-induced resonances and the collective atomic recoil laser," 
Phys. Rev A 59, 585-896 (1999). 

17. H. Y. Ling, H. Pu, L. Baksmaty, and N. P. Bigelow, "Theory of a collective atomic recoil laser," Phys. Rev. A 
63, 053810(2001). 

18. D. Kruse, C. von Cube, C. Zimmermann, and Ph.W. Courteille, 'Observation of lasing mediated by collective 
atomic recoil," Phys. Rev. Lett. 91, 183601 (2003). 

19. M. Gangl and H. Ritsch, "Cold atoms in a high-Q ring cavity," Phys. Rev. A 61, 043405 (2000). 

20. C. von Cube, S. Slama, D. Kruse, C. Zimmermann, and Ph. W. Courteille, G. R. M. Robb, N. Piovella, and and R. 
Bonifacio, "Self-synchronization and dissipation-induced threshold in Collective Atomic Recoil Lasing," Phys. 
Rev. Lett. 93, 083601 (2004) 

21. A.T. Black, J.K. Thompson, and V. Vuletic, "Collective light forces on atoms in resonators," J. Phys. B: At. Mol. 
Opt. Phys. 38, (2005). 

22. P. Domokos and H. Ritsch' "Collective Cooling and Self-Organization of Atoms in a Cavity," Phys. Rev. Lett.89, 
253003 (2002). 

23. J. K. Asboth, P. Domokos, H. Ritsch and A. Vukics, "Self-organization of atoms in a cavity field: Threshold, 
bistability, and scahng laws," Phys. Rev A 72 (5) 053417 (2005). 

24. D. Nagy, J. K. Asboth, P. Domokos, and H. Ritsch, "Self-organization of a laser-driven cold gas in a ring cavity," 
EuroPhys Lett. 74(2), 254 (2006). 

25. G.R.M. Robb, N. Piovella, A. FeiTaro, R. Bonifacio, Ph. W. Courteille, and C. Zimmermann, "Collective atomic 
recoil lasing including friction and diffusion effects," Phys. Rev. A 69, 041403 (R) (2004). 

26. S. Slama, C. von Cube, B. Deh, A. Ludewig, C. Zimmermann and Ph. W. Courteille, "Phase-sensitive detection 
of bragg scattering at ID optical lattices," Phys. Rev. Lett. 94, 193901 (2005). 

27. Th. Elsasser, B. Nagorny and A. Hemmerich, "Optical bistabiUty and collective behavior of atoms trapped in a 
high-Q ring cavity," Phys. Rev A 69, 033403 (2004). 

28. M. Perrin, Z. Ye, and L. M. Narducci, "Microscopic theory of the collective atomic recoil laser in an optical 
resonator: The effects of collisions," Phys. Rev. A 66, 43809 (2002). 

29. C. von Cube, S. Slama, D. Kruse, C. Zimmemiann, Ph. W. Courteille, G. R. M. Robb, N. Piovella and R. 
Bonifacio, "Self-Synchronization and Dissipation-Induced Threshold in Collective Atomic Recoil Lasing," Phys. 
Rev. Lett. 93, 083601 (2004). 



1. Introduction 

Light forces have become an essential tool to cool, trap and manipulate atomic gaseslU. Ulti- 
mately it allowed to reach BEC and thus the lowest temperatures experimentally available||2l. 
Ever since this success there has been an ongoing quest to extend laser cooling to more atomic 
species and also molecules. Here the need for a closed optical cycling transition and the prepa- 
ration of an initial sample of sufficient low temperature and high phase space density to start the 
cooling process has posed severe limitations. In particular for molecules, spontaneous emission 
induced by the laser fields populates a great manifold of electronic and rovibrational levels, 
which hampers the standard laser cooling techniques. 



Recently is has been theoretically predictedO] |4] |5] and experimentally confirmed ||6l|7l[8l 
that using high Q resonators for the cooling light fields, the role of spontaneous emission can be 
strongly suppressed in laser cooling. As one further decisive experiment superradiant enhance- 
ment of cavity cooling was recently demonstrated for Cesium atoms in a confocal cavity 
Although the experiments were performed on atoms, this opens very promising prospects on 
molecules as well. Trapping and cooling of atoms in a far detuned ring cavity field was recently 
also achieved||TOl[II]|221. 

As most sources of gaseous cold molecules are fast beams, the second problem to be dealt 
with is stopping the molecules of such a beam in a very short distance within the vacuum cham- 
ber and trap them. A first success in this direction was achieved by using periodic electrostatic 
fields to stop a molecular beam lfT2l or by creating the molecules already cold by laser ablation 
or directly from an atomic magneto optical trap lfTSll . Here we study the prospects of using a 
ring cavity setup as an alternative way to slow down a gaseous beam. The particles kinetic en- 
ergy is converted into light energy using inelastic coherent backscattering. Such a setup, called 
correlated atomic recoil laser (CARL), was historically proposed in analogy to the Free electron 
laser (FEL) as laser gain medium using atoms lfT4i rT5]| . Many theoretical investigations of these 
phenomenon, sometimes also called recoil induced resonance (RIR), produced partly contro- 
versial results and interpretations lfT6ll28l . For quite some time experimental efforts to observe 
clear signatures of CARL were only partly successful fTTl. 

Only recently, a first unambiguous observation of the (inverse-) CARL mechanism to accel- 
erate a cold gas by coherent light scattering in a cavity was achieved in Tubingen using Cesium 
atoms ifTSi 26 il. These results are in very good agreement with predictions from microscopic dy- 
namic model s fT9l as well as collective macroscopic descriptions of atoms in ring cavities ll20l . 
Even stronger stopping forces were observed and extrapolated by Vuletic in a standing wave 
cavity setup 191 1211 . where a special bistable atomic selforganization process superradiantly en- 
hances the light forces If3l l23]| . In all of the experiments the particles were atoms with a closed 
optical cycle, which were cold and relatively slow from the start. Nevertheless these results 
show that the method has great promise. 

Here we investigate the experimental conditions necessary to stop and trap a fast molecular 
beam by help of a single side pumped ring cavity geometry. The idea relates closely to the 
original idea of CARL, but we are interested in stopping and cooling particles and not in laser 
gain. In contrast to most previous work we consider the limit of large detuning, where inversion 
and spontaneous emission from the particles plays no role. This enables one to adiabatically 
eliminate the excited states from the dynamics, which allows us to apply the results to any 
point like objects with suitable polarizability and in particular cold molecules. 

In this limit the equations simplify dramatically but, of course, the atom field interaction 
strength per particle gets very small and only cooperative light scattering by a large number of 
scatters can induce a sufficiently strong force. In the quasi ID geometry along the cavity axis 
this requires a periodic particle arrangement, so that they act like a Bragg mirror coupling the 
light fields within the resonator For favorable parameters the interference of the injected and 
backscattered light creates an optical potential suitable to stabilize such a regular distribution 
or even induce it via selforganization from an initial flat distribution. 

Note that momentum conservation assures that the total force on the cloud is directly related 
to the backscattered field intensity amounting to two photon momenta per photon decaying 
from the backward mode. Hence the average force on the gas will always point to the same 
direction and this acceleration prevents the system from reaching a steady state in the long run. 
Nevertheless at an intermediate timescale one can expect the atoms to form a transient spatial 
equilibrium distribution with respect to their accelerated center of mass. 

In the following we try to answer two closely related questions in this respect: ( 1 ) Under what 



conditions and circumstances is a homogeneous particle distribution within the resonator unsta- 
ble against small density fluctuations so that the atoms will start to selforganize? (2) How does 
the self consistent atom field distribution look like in the regime where self ordering happens 
and what is the corresponding field intensity and average force on the cloud? As a third impor- 
tant question we address the scaling of the system with particle number, volume and required 
pump intensity to see whether the effect persists on a macroscopic scale. 

Methodically we will adapt an idea already successfully applied to selforganization of atoms 
in standing wave cavity fields. For this one iteratively determines the self consistent steady state 
of the cavity field with particles ordered in a way just to create this very same field. The iteration 
is done assuming fast atomic thermalization and field equilibration but only slow change of the 
average atomic temperature. In order to adapt the method to particles in a ring cavity we have 
to change to the accelerated rest frame of the atomic clouds center of mass and additionally 
assume that the equilibration is faster than the center of mass acceleration. This allows us to 
find an accelerating position distribution with a largely time independent shape. Investigating 
the stability and self consistence of this distribution will then give us limiting conditions for 
which the proposed method should work. In order to check the predictions we then will directly 
simulate the microscopic equations of motions for finite ensembles of particles and the fields. 



2. Model 



As model system let us consider a large number of linearly polarizable point particles in- 
teracting with two degenerate counter propagating linearly polarized plane wave ring cavity 
modes. 

The electric field inside the cavity is then 



'E{x,t) = t'^+\x,t)e-''"'" + t^-\x,t)e"^'" , 



E'^+\x,t) = £ a+{t)e''' + a^{t)e-'''' 



(1) 
(2) 

(3) 



where the electric field per photon is e = ^Jhaip/leoV and a±{t) gives the field amplitude. 

The pump and cavity frequency are assumed to be sufficiently detuned from any optical 
resonance so that the particles are only weakly excited and the polarization can be adiabatically 
eliminated. Thus they experience a Stark shift linear in the local intensity and act like a dynamic 
refractive index in the field. Let the parameter Uq denote the light shift per photon in the field. 
From energy conservation one sees that Uq then also has to give the frequency shift of the cavity 
mode resonances per atom. In terms of the two level linewidth F, detuning A„ and atom-field 
coupling g the effective mode frequency shift reads: 



F2+A2^ 



(4) 



In a typical setup one sends the molecular beam at average velocity vq into the ringcavity 
nearly antiparallel to the driven circulating field. The molecules might be moving quite fast (up 
to "100 m/s) without some pre-slowing so that the ratio of the Doppler shift to the linewidth 
kv/ can be larger than unity but not too large. At much larger velocities, the pump and or the 
backscattered fields are far of resonance, which requires stronger and stronger pump amplitudes 
to generate a noticeable force. 



In the following we will shift our reference frame to coincide with the particles center of 
mass motion. Hence, while a particle at rest in the cavity would see a single cavity detuning 
Ac, the molecular beam experiences two different Doppler shifted cavity fields with detunings 
A± = Ac ± kv. Thus the equation of motions for the field mode amplitudes are lfT9l : 

^fl±(f ) = [/A± - iNUo - k] fl±(f) - iNUo ^eT2to^ ^^(f) + (5) 

Here denotes the particle number and rj^ gives the cavity pump amplitudes. We see that 
the particles density distribution couples the two counter propagating modes via the spatial 
average R± — (e'^^'*^''), which is called particle bunching parameter and one has ~ A 
uniform distribution gives R± « 0, while an optimally ordered phase gives \R± | ^ 1 . 

As second part in our model we include the force of the light field on the particles. Again 
assuming the large detuning limit this force is dominated by the dipole force lfT9| . This force can 
be derived from the effective optical potential, which is determined by the local field intensity 
simply by: 

e^V{x) = Uo\E{x)\\ (6) 

For Uo < 0, i.e. Aa < 0, this attracts atoms to maxima of the cavity field. Using the definition 
of the field in terms of the amphtudes a±, E^{x) ~ e(fl+e''"^ + fl-e^'^'^) we get: 



= + |fl_P+2|fl+fll|cos(2ytx + a) 



(7) 
(8) 



with a being the argument of a+a*_. Note that for calculating the force the constant terms can 
be omitted from the potential and the average force then simply reads: 

f{x) = -UqV \E{x)\^ = -2TikUQ\a+a*_\sm{2kx + a). 

This force will then act on the atomic distribution and in general one gets a rather com- 
plex coupled atom field dynamics [ 22 1. Note that neglecting spontaneous emission (no radiation 
pressure force) we get momentum conservation of the cavity field modes and gas. This has the 
consequence that for single side pumping, for every photon in the counter propagating mode 
one has transferred a momentum of 2hk to the cloud. Hence in steady state the average force 
is directly proportional to the output intensity of the counter propagating mode and we will 
always have acceleration in the pump direction. 



3. Coupled particle field dynamics 

3.1. Quasistationary field 

Lets us now get some first qualitative insight in the ongoing dynamics for a thermal cloud in a 
single side pumped ring cavity (77^ = 77 , 77 _ = 0). While for an empty cavity the field intensity 
and thus the optical potential is constant along the axis, the particles will coherently backscatter 
some light and create a periodic intensity modulation. In a typical experimental setup this field 
redistribution will happen much faster than the typical time on which the particle distribution 
evolves. Hence the field will almost instantaneously reflect the momentary particle distribution. 
As the only relevant quantity of the particle distribution entering the field equations is the 
parameter R = (e^''^^^, the field amplitudes will relax towards: 



iA+-NUo+iK){A^^NUo+iK)-N^U§ |/?_| 2 ' ^ ^ 

^. NUqR^ 

'^(A+-A^f/o+/K-)(A_-A^f/o+'K-)-A^2t/o 

To get some more insight in the physical meaning of R let us consider the prototype case 
of a periodically modulated gas density with symmetric maxima at x = xq + nX /2. R than can 
be written as R ~ |/?_ | e^'fao^ where the modulus |/?_ | measures the modulation depth and the 
location of the maximum density determines the phase of R . Inserting these expressions into 
the field equations, we then get; 

* 1^,2 (A_-jV^o+/>c)jV^o|/;_|g-'fao 

fl+fl_ = |T7| (11) 

I (A+ -NUq+Ik) (A_ -NUo+iK)-N'-U^ \R-\^\ 

Here we see that the better the atomic localization | the more backscattering and stronger 
modulation of the intracavity intensity we get. As a second important piece of information 
this result allows us now to determine the spatial shift a of the cavity field intensity maxima 
with respect to xq by help of Eq]?] Since Uq = ~ \Uq\ = |?7o|e"^ we define the relative shift 
5.x = ^ /(2A:) of the field intensity maxima with respect to xq by 



A_ -NUo + iK = ^(A_ -NUof + K^e''^"'-'"' . (12) 

This gives a± = (p — 2kxo = —2k{xo — (p /(2fc)) = —2k{xo — 8x), where the value of (j) is deter- 
mined by (fTTT i: 

-A^+NUo . ^ -K 

cos 9 = — z, sm0 



v/(A_-A^C/o)2 + ^2' yiAT^^MToF+i? 

As the maxima of the light field intensity determine the minima of the optical potential Eq. 
|6] (stable equilibrium points), we see immediately that field minima induced by the density dis- 
tribution of the gas do not coincide with the peak density positions xq but are shifted by 5v. 
Hence we obviously cannot find a self consistent time independent configuration of particles 
and field as for the standing wave case i23| . It is very interesting to note though, that the magni- 
tude of the shift only depends on the total particle number via NUq and cavity parameters, but 
not on the pump strength 77 or the precise form of the atomic distribution beyond the average 
(e2too^. 

Obviously the above calculated shift of the potential implies an average force on the cloud 
and an acceleration of the particles center of mass position. This result can also be derived 
from basic momentum conservation arguments as follows: since all photons in the backscat- 
tered mode have lead to a deposition of 2hk of momentum in the cloud, the gas is accelerated, 
whenever a backscattered field is present. As a consequence of this acceleration, R will be 
time dependent in general. Hence in contrast to atomic selforganization in a standing wave 
cavity field||24|, we cannot expect the system to reach a time independent self consistent atom- 
field steady state. Nevertheless, as we will see below, quasistationary self consistent atom field 
distributions can still form, if this acceleration is not too fast. 

3.2. Quasistationary accelerated atomic distribution 

As argued above the coupled atom field dynamics produces an acceleration of the particles cen- 
ter of mass. As long as this acceleration is not too fast the corresponding periodic field intensity 



distribution will almost instantaneously follow the particles center of mass motion. Hence the 
particles see only a slowly time varying optical potential relative to their center of mass position. 
This should allow the particle distribution to adjust to this accelerating potential. In order to get 
some quantitative insight here we will now transform our equations to a suitable accelerated 
reference. Primarily this introduces an inertial force (acceleration): F,„ = mg into the atomic 
equations of motion, which can be accounted by adding a term V,,, = —mgx to the effective 
atomic potential. Obviously a suitable choice of g requires to shift the periodic local potential 
minima to coincide with the peaks of the atomic position distribution. This allows subsequently 
to find an approximate stationary atomic distribution with respect to the accelerated potential. 

Of course the transformation to this frame is reflected in the field equations, where the de- 
tunings are now time-dependent quantities induced by the time dependent Doppler shifts of the 
cavity modes©. Nevertheless, as long as the acceleration is not too fast on the timescale of the 
field evolution if, the fields will still adiabatically follow changes of atomic velocity and dis- 
tribution. Lets consider this in a bit more detail now. Omitting the constants parts, the potential 
in the accelerated frame reads: 

V{x) = - \Uo\ \a+a*_\cos{2k[x-{xQ--^)])~mgx, (13) 
' ' 2k 

where we choose ^ in a way, that the peaks of the atomic density distribution XQ + nX/2 coincide 

with local potential minima, i.e. we want (d /dx)V{x) = at x = xq + nX /2. This requires 

~2k\UQ\\a+a*_\sm(p +mg = 0. (14) 

Substituting the steady state expressions for sin (j) and for a+a*l, we then get 



2k a^|;7o|^|/?_|k- 



g^^\n^\— " '•' y. (15) 

\{A+-NUQ+iK)iA--NUQ+iK)-N^U^\R-\^\ 

The resulting acceleration g thus is negative and depends on the magnitude of the atomic 
locahzation parameter R . This is consistent with our expectations, since the mode e^'*^^ is 
pumped, and the scattering converts photons travelUng the negative x-direction to photons in 
the positive direction and we have 

k 

mg = -4K\a*_a^\—. (16) 

Besides providing the basis for the choice of an accelerated frame, Eqs{T5][l6]are quite inter- 
esting and useful results on their own. While for small the acceleration linearly grows with 

showing the superradiant A^^ enhancement of backscattering, for large enough it decays 
like 1 /N. Similarly for large K (bad cavity) and a nearly flat distribution R ^ acceleration 
vanishes. In these limits we can expect the distribution to have enough time to equilibrate to a 
quasistationary distribution in the accelerated frame. 

Note that the acceleration for not too large A^ depends on \Uq \ and thus on the inverse square 
of the particle-field detuning as it is the case for near resonant radiation pressure and sponta- 
neous emission. Hence at first sight this does not look favorable for molecules. However, the 
coflective enhancement factor A^ in the force will not show up in spontaneous emission, which 
allows to strongly increase the force per spontaneously scattered photon in the resonators. 

Having now found a proper frame as a next step, we will calculate a corresponding atomic 
distribution. For this we adapt the ideas developed for the standing wave cavity case ll23l and 
simply use a thermodynamic equilibrium canonical density distribution for the particles. Hence 
the particle distribution is approximated by 

pix) = yxp[-V{x)/kBT], (17) 



where the normalization is given by the partition function 



J exp[~V{x)/kBT]dx. (18) 

Note that the effective potential now contains the acceleration terms mgx and this procedure 
only makes sense as long as V{x) has periodic local minima of sufficient depth. This is definitely 
be the case for small g as it occurs for an initial nearly homogeneous particle distribution, where 
R is close to zero. Hence we can at least use this formalism to study the initial phase and 
threshold conditions of the selforganization process. 

3.3. Selforganization threshold 

As we have seen above the cavity field and inertial force from the acceleration g determine 
an effective combined potential V{x) for the accelerated particles for which we get a corre- 
sponding stationary distribution from Eq[T7] This modified density distribution corresponds to 
a new bunching parameter/? and a new optical potential. By subsequent application of Eqs.fTTl 
[13] this procedure can be iterated and in many cases a self consistent solution is obtained for 
p(x) and V{x) after a few iterations. Note that in order to get a meaningful result, we have to 
restrict ourselves to a single potential well (local potential minimum) to determine the density 
distribution p{x) and assume that all other wells to have a similar number of particles. This 
can be physically motivated from the fact that redistribution between different wells will occur 
on a much longer time scale as internal thermalization. Note the we do not account for direct 
interparticle interactions here, which would occur for very high local densities. 

There are now several decisive aspects that can be learned about the physics of our problem 
from this mathematical iteration procedure. We will start here with the most simple but central 
question, whether a homogeneous initial particle distribution at rest within the resonator is un- 
stable against small perturbations. This amounts to a CARL type setup, where particles, which 
are originally evenly distributed are collectively accelerated by superradiant scattering ll25l . In 
the following we will first determine the minimal pumping threshold required to initiate such 
superradiant backscattering via self ordering. In a second step we will then check, whether a 
quasistationary accelerated ordered phase can develop out of these initial fluctuations. 

Starting from an initial flat particle distribution, we obviously have R_ ~ and the atoms 
generate no coherent backscattered field. Hence the corresponding average acceleration is zero 
at the beginning and we can start in laboratory rest frame for the dynamics. Superradiant back- 
scattering and collective forces in this case can only build up, if this solution is unstable against 
small density fluctuations with period X/2. Hence we add an infinitesimal perturbation to our 
uniform distribution p (x) = j in the form: 

p(o)(x) = i(l+er(x)). (19) 
r{x) in general could be any periodic function 

= £ iA,„cos{mkx)+B,„smimkx)), £ aI+bI = 1 (20) 

m= 1 m= 1 

and see whether this perturbation is amplified during an iteration cycle. Note that within the 
ring cavity geometry the only perturbation that contributes to R is that for m — 2, i.e. we 
can limit ourselves to r{x) — A2Cos{2k{x — XQ))+B2sm{2k{x — xq)) as a perturbation. Due to 
translational invariance of the problem we can also set B2 ~ 0. Hence the perturbation simply 
leads to R — eA2/2 for the mode coupling parameter. 



As a first step in the iteration we will now use this R to determine the corresponding fields 
and the optical potential. As shown above in Eqll3]the optical potential is shifted with respect 
to XQ by an amount independent of A2. In the second step we now choose the acceleration ^ in a 
way described above to pin the effective potential minimum to xq and apply Eqs . [T3] and [TT] to 
find the corresponding thermal density distribution for this potential. Of course this distribution 
now shows a modulation induced by the optical potential which we have to compare with the 
magnitude of our initial perturbation. 

As e is small we do this analytically by expanding p (x) in Eq(T7]to lowest order in e to find: 



V{x) 



-(1 



eA2 



)/A 



NUi\ri\ 



ksT \{A^-NUo-iK){{A+-NUo) 



(21) 



■[cos{2k{x-XQ)] + 0{e^) (22) 



The central question now is, whether the density perturbation has been enlarged or dimin- 
ished in this iteration cycle. For a flat distribution to be stable, the first order correction term in 
e must be less than A2. Thus stability requires: 



^thresh ^ 



\ 



kBT^/iA^-NUof + K^ ((A+ -NUo)^ + K^) 



(23) 



This formula is one of the central results of this work, which will discuss in more detail at 
different examples in the following. It allows to clearly determine the experimental conditions 
for the cavity, initial particle phase space density and laser power required for such setup. We 
will later also substantiate the key properties of this formula by direct multi particle simulations. 

4. Results and discussion 

4.1. Threshold condition and acceleration of a gas at rest 

To investigate the underlying physics in more detail let us first consider a cloud at rest v = 0, 
i.e A+ = A_ = Ac. In this case the threshold condition reads: 



V thresh ^ 
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(24) 



Obviously a flat distribution always gets unstable for sufficient pump power. One needs less 
pump intensity the more particles one has and the stronger they are coupled to cavity mode 
(bigger modulus of Uq). The required pump intensity scales linearly with the initial temperature, 
which gives some important restriction for useful molecular sources in particular Of course 
the validity of the threshold formula has its limitations as the weak saturation and two level 
approximation for the particles eventually will break down at certain pump intensities. 

Note that this stability threshold for a flat distribution attains a minimum value, when A^ — 
NUq, where it simplifies to 
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(25) 



Of course A^- = NUq is just the resonance condition for the cavity mode shifted by the refrac- 
tive index of the particles and hence gives a maximum of the intracavity photon number. As we 





Fig. 2. Same as Fig.Qfor A' = 200, Uq = -1/115, K = 1. 



will see later from Eq l3.1l this is a rather tricky parameter choice. In this case the corresponding 
shift of the optical potential minima with respect to the atomic center of mass just corresponds 
to 7r/4. Here the particles sit on the position of maximum slope and the optical dipole force has 
maximum as well. However, as the curvature of the potential is zero at this point the effective 
potential in the accelerated frame has no local minimum and one cannot expect the formation 
of a self consistent atomic distribution of the form of EqlTTl 

In the following we will investigate the validity of this threshold condition numerically by 
iterating equations 1131 1 7l for a selfconsistent atom-field configuration. We see that the iteration 
results agree nicely with the analytic condition of Eq|25] As shown in Figs. II 121 l3]andl2]we find 
a numerical instability in agreement with Eq|25]with the minimum threshold at = NUq. The 
onset of self ordering contour is readily seen to closely follow the red dashed line representing 
the analytic threshold condition Eq|24| = 200, C/q = -1/115, «•=!. 



Contour of max density from itreation simulation and threshold condition 
forN=1000 and Uq=-0.002k 
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Fig. 3. Same as Fig.EforA' = 1000,(70 = -1/375, K=l. 



Simulated kinetic energy gain for 200 particles 




Fig. 4. Contour lines of kinetic energy gain of N=200 particles with average initial velocity 
zero and average kinetic energy hK and Uq = — 0.015 ic within a fixed time interval Kt = 60. 

Let us now further substantiate these predictions for some special parameter choices by a 
real space dynamical solution of the equation of motion of N particles and the field modes. 
For this we use a semiclassical dynamical model as developed in ref. ISl, where the atoms are 
modeled by polarizable point particles and the fields are amplitudes are represented by coherent 
states. This model has already been successfully applied to simulate cavity enhanced laser 
cooling and atomic selforganization in a standing wave cavity geometry ll23l . Starting from a flat 
thermal distribution atomic distribution and with zero field ampUtude, we then simply monitor 
the coupled atom-field evolution for a set of typical parameters. Due to the large number of 
variables we have to restrict ourselves to monitor only some characteristic observable like the 
evolution of the total kinetic energy, the magnitude of the bunching or the time evolution of the 
center of mass of the cloud. This is shown in Fig|4]where we compare contour lines the kinetic 
energy gain of the gas after a time of Kt = 60(check) with the analytic threshold formula from 
above. Indeed selfordering and significant acceleration of the particles is found only above the 
analytic threshold line. 
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Fig. 5. Self consistent atomic density distribution below and above threshold for A' = 2000 
atoms for Uq = — 0.001 ic and Sc = —2k. The analytic threshold result for this set of pa- 
rameters is rjiiji-g^ij = 31.6. 

4.2. Quasistationary accelerated particle distribution 

As we now found a condition for the breakdown of stability of a homogeneous distribution, 
the next big question now is what happens in this instable regime? What kind of dynamics can 
we expect in general and under what conditions can we expect selforganization and an ordered 
accelerated particle distribution? 

Obviously the backscattered field amplitude gets nonzero and the center of mass of the cloud 
gets accelerated. As a first guess for the further evolution we will now simply iterate the above 
loop further and see under which conditions a self consistent accelerated atom field distribution 
could possibly form. As we have no reliable time scales for the atomic thermalization process 
available, we will simply assume that it is faster than the time scale of significant acceleration 
of the cloud so that we keep the average velocity constant in these iterations. 

Experimentally and in numerical simulations it has been found that in particular the regime 
NUq > K has the biggest potential for a self ordered phase and interesting dynamics, so we will 
primarily focus on this regime 1231 . 

In general for a wide range of parameters a self consistent accelerated atom field distribution 
can be found at a sufficient pump strength. A plot of the density distribution below, at, and 
above pump threshold is shown in Fig |5] Actually the iterations show, that the onset of self 
ordering depends non trivially on all the parameters and clear indications can already be seen 
even after only about five iterations. Typical numerical results can be visualized by plotting 
the peak atomic density versus pump strength for 5, 10 and 100 iterations as in Fig|6l There 
are no significant differences in convergence between 50 and 100 iterations. This peak density 
is a good indicator of atomic ordering, which due to translational invariance may form at any 
location between the mirrors. Alternatively the self consistent backscattered intensity turns out 
to be a good indicator for the onset of self ordering as a function of pump strength. From this 
quantity we can also determine the average force on the cloud. 

The self ordering threshold can be located in Fig|6]or from the contour plots such as Fig|2] 
or [5] via the transition from the flat homogeneously distributed background to the exponential 
increase of local particle density. The transition is easy to see by eye. Numerically we have to 
select some value of particle density larger than the homogenous density and above the "noise 
background" as an indicator of selforganization. In this way we can generate contour lines of 



Peak density distribution and acceleration factor g versus pump ampiitude. 
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Fig. 6. Maximum atomic density as a function of pump amplitude for three different itera- 
tion numbers of 5 for dotted line, 10 for dashed line, and 100 for solid line. 



the maximum density of the self consistent particle density. This allows a good comparison with 
the analytic form for the pumping threshold as a function of other parameters such as detuning. 
A plot of the density distribution below, at, and above threshold is shown in Fig ID The key 
feature to note in Figs [2] and [3] is that in accordance with the analytical results derived above, 
the minimum threshold which occurs, when = NUq, with a symmetric increase around this 
minimum. In seems that for a large range of parameters the instability of the flat distribution 
is connected with the existence of an ordered accelerated phase which is particularly stable for 
A, < NUo. 

A decrease of cavity damping leads to a sharper resonance of backscattered light at certain 
values of cavity detuning. We can estimate the value of cavity detuning which leads to the 
maximum backscattering (and average force) while holding other parameters fixed by taking 
the derivative of |fl_ 

A, F,na. = NUo ± ^Jn^U^\R-\^-k^ (26) 

Recall the bunching parameter approaches unity here. This result is only useful for small 
damping, k < NUq, where there is a distinct, sharp maxima of |fl- 1^ in the numerical results at 
Ac = and A^ = 2NUo (with very little backscattering at A^ = NUq). Thus as seen in Fig. [7] 
there are effectively these two resonant frequency of the cavity, with a larger Q leading to more 
backscattered light, but over a narrower A^ frequency interval. 

Alternatively, for larger values of damping, k > NUq, the numerical results for backscattering 
show a maximum value at Ac = NUo which decreases slowly away from this value of detuning. 
The limit on the overall range of cavity detuning which still permits self ordering of the particles 
is set by the pump. The threshold pump value given in Eqs [25] and l24l should be exceeded to 
observe these backscattering dependencies on detuning. 

To get some quantitative comparison with the selfconsistent atomic distributions calculated 
by the above iteration procedure and direct numerical simulations, we now also depict a contour 
plot the maximum of the atomic density distribution of a cloud of N ~ 1000 atoms NUo = 
—2k after a given time interval At = xx/k . We see that clear signs of selforganization appear 
when we cross the theoretically found threshold condition depicted as the line with crosses as 
a function of detuning and pump strength. We also see that crossing the threshold and ordering 
is accompanied by acceleration as expected. 



Contours of back scattered intensity from iteration simulation 
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Fig. 7. Intensity of backscattered field mode as a function of detuning 8c and pump ampli- 
tude T] for A' = 1000 atoms Uq = -0.002k: 




Contour of density from iteration simulation and ttirestioid condition 
for N=3000, U^=-0.0027, and l<v=5. 
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Fig. 9. Contour lines of maximum of the self consistent atomic density for N = 3000 atoms 
for (/() = —0.0017k" and k\' = 5k. For comparison we also show the analytical expression 
for the self organization threshold amplitude multiplied by 3 (dotted line). 

4.3. Threshold conditions and deceleration prospects for a moving thermal cloud 

Having investigated the dynamics of a gas at rest in some detail, we now come back to one of 
our original central goals and check the possibility of trapping and cooling molecules from a 
beam by superradiant scattering. In a typical setup one would send a molecular beam through 
the cavity in a small angle almost parallel to the ring cavity axis and inject a strong counter 
propagating coherent field. In such a setup the molecules are initially moving quite fast (up 
to ~100 m/s) with constant spatial density. The central question now is for which phase space 
densities, velocities and pump amplitudes superradiant scattering will appear and can be used 
slow down and trap the molecules. 

In general using the rest frame of the beam the equations are the same as above except for 
different detunings A± = ± kv. for both cavity modes. Hence one can expect a similar thresh- 
old behavior. Naturally for a large value of the ratio of kv/K the pump and the backscattered 
field cannot be on resonance simultaneously. This will increase the threshold or even make self 
ordering not be possible. In principle on has the two possibilities of resonant pumping but a 
non resonant backscattered field or off resonant pumping to get resonant enhanced of the back- 
scattered field. The threshold for self ordering is first calculated assuming the cloud velocity is 
constant throughout the iteration. We can compare this result to the threshold obtained using 
the accelerated frame where the velocity is updated throughout the iteration. 

The most notable feature of the numerical results as seen in Fig.|9]is the asymmetry in the 
pumping threshold with detuning as compared to the case with particles at rest. 

Additionally, the analytic stability analysis leads to a first order stability equation for optical 
pumping threshold given by Eq.l23] 

For NUq > K there is once again excellent agreement between the numerical and analytic 
results for pumping threshold as shown in Fig. |9]with only about a two percent difference in 
threshold values. Trapping molecular beams with relatively large velocities is possible with a 
large pump and large cavity decay rate. In this region where NUq < K, the analytic and numer- 
ical results for threshold are less in agreement, although the results have the same functional 
form as kv increases. The backscattered intensity as a function of pumping versus velocity of 
the moving cloud is shown in Fig.|7] 



Bunching parameter |R| for N=3000 and p„=2000hk 
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Fig. 10. (a) Kinetic energy loss per particle within a time interval of Kt = 60 as a function 
of detuning and pump amplitude Tj for = 3000 atoms counterpropagating the cavity 
axis with initial average momentum of po = 20007; ff and temperature kJQ = TiK and (b) 
average bunching parameter during this time for the same parameters. 

4.4. Numerical multiparticle simulations of cold beam stopping 

As a final check and demonstration of the efficiency of the setup we will now present some 
results on the possibility to stop a fast beam in such a setup again using real space simulations. 
In this case we start with beam of N=3000 particles in a thermal distribution of temperature 
kT = hK but with additional average momentum of p- = — 2000M per particle in the negative 
z-direction counter propagating the pump. As a measure of the efficiency of the scheme we 
first plot the amount of average kinetic energy extracted from the atoms by the fields in a time 
interval Kt = 60 as a function of pump strength 77 and pump detuning . 

Obviously we get efficient deceleration of the cloud when either the pump mode directly or 
the reflected light field is in resonance with the effective cavity eigenfrequencies ifTol . When we 
plot the time averaged bunching parameter we see that efficient energy transfer is associated 
with atomic selfordering as analytically predicted. As the atomic bunching is more efficient for 
negative detunings we get a more efficient energy extraction in this limit. 

These results compare qualitatively quite well with the analytic prediction of Eq|23]as shown 
in Fig |4.4l where we compare a contour plot of bunching parameter with the analytic threshold 
formula for the selforganization threshold multiplied by 3 (line with crosses). 

5. Conclusions 

We studied the conditions to stop and trap a fast molecular beam by help of a single side 
pumped high Q ring cavity geometry in a CARL setup. The analytical results for kinetic energy 
extraction from a self consistent stationary solution of the coupled particle field equations were 
checked by numerical simulations. In the limit of large detuning inversion is negligible and 
coherent light scattering dominates. As the atom field interaction strength per particle is small in 
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Fig. 11. Contour lines of time averaged bunching parameter \R\ found by numerical simu- 
lation as function of detuning Sc and pump amplitude eta for parameters above. Again for 
reference we show the analytical result for the selforganization threshold (x3) 

this regime we need collective enhancement of the Raleigh scattering blue shifted with respect 
to the excitation light to do the job. This can be achieved in a regime where the particles arrange 
in a periodic density lattice and form a moving Bragg mirror. 

Finding the effective optical potential in an accelerated frame and applying a linear stability 
analysis allows to derive the analytic threshold for the laser pump power, above which an initial 
homogenous particle distribution gets unstable against small density fluctuations. We see that 
the available phase space density at the beginning determines the required threshold laser inten- 
sity. Above this threshold a suitable choice of detunings causes the interference of the injected 
and collectively backscattered light fields to establish a periodic particle density acting like a 
moving Bragg mirror, which yields a time dependent coupling of the counter propagating fields 
within the resonator. By using a mean-field analysis the self consistent atom-field distribution 
is found by an iteratively which allows to quantitatively predict the maximal acceleration and 
localization of the cloud. 

In a second step we checked these predictions by simulating the microscopic equations of 
motions of the particles for finite ensembles. These simulations show surprisingly good quali- 
tative and quantitative agreement with the analytic threshold result and acceleration predictions 
within the regime, where the analytical model predicts selforganization. Outside this regime, 
where the effective potential in the accelerated frame shows no periodic local minima the sim- 
ulations differ from the analytic predictions. This can be traced to cavity induced heating pro- 
cesses in this regime as well as remnant radiation pressure not accounted for in the analytic 
model. In any case the analytically found threshold formula for the selforganization process 
provides a solid criterion for the prospects of a corresponding experimental setup. 

As a bottom line trapping and stopping of a fast molecular beam in a ring cavity seems exper- 
imentally feasible with current technology, provided one can find a strong enough pump laser 



at a frequency where the particles are polarizable but sufficiently transparent. The higher the 
initial phase space density of the particles, the faster this process will happen and the stronger 
are the deceleration forces and the suppression of spontaneous emission. In favorable cases this 
could be combined with cavity cooling in the particles rest frame as well|3i|. 

This work could be extended by continuing the work of Asboth who has used the mean 
field analysis to find self-organization thresholds in standing wave cavities, but now analyzed 
with transverse mode degeneracy. In particular, Vuletic's group has found the degeneracy of a 
confocal resonator significantly enhances the cooling efficiency. It would be useful to examine 
such threshold behavior in this highly degenerate regime. 
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